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When a surface wave interacts with a vertical vortex in shallow water the 
latter induces a dislocation in the incident wavefronts that is analogous to 
what happens in the Aharonov-Bohm effect for the scattering of electrons by 
a confined magnetic field. In addition to this global similarity between these 
two physical systems there is scattering. This paper reports a detailed calcu- 
C3 [ lation of this scattering, which is quantitatively different from the electronic 

case in that a surface wave penetrates the inside of a vortex while electrons 
do not penetrate a solenoid. This difference, together with an additional dif- 
ference in the equations that govern both physical systems lead to a quite 
different scattering in the case of surface waves, whose main characteristic is 
a strong asymmetry in the scattering cross section. The assumptions and ap- 
proximations under which these effects happen are carefully considered, and 
their applicability to the case of scattering of acoustic waves by vorticity is 
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noted. 

03.40.Kf, 47.35.+i, 47.10.+g 



Typeset using REVTgX 



I. INTRODUCTION 



In a remarkable paper, Berry et. al. [Q clarified the way in which a curl-free magnetic 
vector potential modifies the wavefront structure of an electronic wavefunction that obeys the 
non relativistic Schrodinger equation. They concluded that for electrons travelling outside 
an infinitely long cylinder enclosing a magnetic field, the wavefronts outside the cylinder 
would be dislocated by an amount proportional to the amount of magnetic flux within 
the cylinder. Reasoning by analogy, they also concluded that such dislocated wavefronts 
should occur for surface water waves when they encounter a vortex. A simple experiment 
conclusively demonstrated this effect IJ. 

In the case of the electronic wavefunction interacting with a confined magnetic field 
(and its unconfined vector potential) Berry et. al. |]] also computed the complete solution 
to the Schrodinger equation that, in addition to the dislocated wave, includes a scattered 
wave. Trying to do this in the case of the water waves is however more difficult because 
the analogy between de Broglie waves and water waves breaks down when pushing it into a 
quantitative statement. There are two essential differences: The first is that for an electron 
the appropriate boundary condition is that the wave function vanishes at the surface of the 
cylinder; in the case of water waves, the waves of course penetrate inside the vortex and it 
becomes necessary to solve the appropriate equations not only outside the vortex but also 
inside, and match them with continuity conditions. The second is that the wave equations 
that govern both phenomena, although similar, differ in quantitative details. This paper 
addresses both these issues. 

The scattering of surface waves by vertical vorticity in shallow water was discussed by 
Cerda and Lund || and by Umeki and Lund |J who discovered that a vortex may support 
spiral wave solutions. Fabrikant and Raevsky [|]] have studied the case of a fluid of arbitrary 
depth in a Born approximation. The interaction between surface waves and vertical vorticity 
is in many respects similar to that of acoustic waves and vorticity, a topic that has been 
much studied over the years and that recently has been the subject of particular interest 
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due to the possibility of using acoustic waves as a nonintrusive probe of vortical flows, both 
laminar and turbulent ||, much in the same way that X rays and neutrons are used to probe 
condensed matter structures, both ordered and disordered. Most treatments, however, rely 
on a Born approximation || whose validity breaks down when a surface wave interacts with 
a vertical vortex with nonvanishing circulation, leading to a long range velocity field that 
decays like 1/r where r is the distance to the vortex core. This is precisely the case that is 
studied in the present paper. 

This paper is organized as follows: Section 2 derives the equations that describe the 
scattering of a surface wave by vertical vorticity in shallow water. We pay particular at- 
tention to the assumptions needed to derive those equations. Section 3 has a reminder on 
the Aharonov-Bohm effect as relevant for the present discussion. Section 4 presents the 
computation of the solution to the equations derived in Section 2. Section 5 presents several 
illustrative examples and Section 6 has concluding remarks. Technical details are contained 
in two Appendices. A subsequent paper |7j studies the first order corrections to the shallow 
water approximation. 

II. SHALLOW WATER WAVES IN INTERACTION WITH A VERTICAL 

VORTEX 

We consider the problem of the interaction of shallow water surface waves in an inviscid 
incompressible fluid of uniform depth h with a stationary vertical vortex. The coordinates 
are (x, y) = x in the horizontal direction and z in the vertical direction. The velocity 
and the surface displacement are denoted by v(x, z, t) = (v±(x, z, t),w(x, z, t)) and i](x, t), 
respectively. 

The equation of motion is || 

d t v + (v ■ V> = -p' 1 Vp - g, (2.1) 

where p is the (uniform) density of the fluid, p is the pressure, g is the gravitational acceler- 
ation, d t = d/dt and V the three-dimensional gradient. The boundary conditions on w are 
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w = at the fluid's bottom (z = 0), and 

w = d t r] + v± ■ V±r) (z = h + r]{x,t)), (2.2) 

where Vj_ is a horizontal gradient, at the surface. We consider a free surface and neglect 
surface tension, which is consistent with the shallow water approximation. 

In shallow water the length scale for spatial variations is much bigger than the fluid depth 
h. Consequently, the continuity equation V ■ v = together with the boundary condition 
at the bottom imply 

w(x, zt)\ z=h+v = -Vl • v ± \ z=0 {h + r,) (2.3) 

to leading order in h/L, where L is the length scale for space variations. Inserting ( |2.3| ) 
into the kinematic boundary condition ( |2.2| ) and using that, to leading order in the shallow 
water approximation 

v±\ z=0 = v±\ z=h+v 

leads to 

d t T] + hV x -v ± + V±- (r]v ± ) = 0, (2.4) 

assuming surface deformations small compared to depth (77 <C h). 

Neglecting vertical accelerations with respect to g, the z-component of fl2.1|) and conti- 
nuity of the pressure at the free surface yield 

p(x, z, t) = pg(h + r](x, t)-z)+ p a , (2.5) 

where p a is the atmospheric pressure. Substitution into the cc-component of ( j2.1|) gives, 
again in the shallow water approximation 

d t v± + (v± ■ S7±)v ± = -gV±r]. (2.6) 

We will consider surface waves with particle velocity u(x, t) and surface deformation 
r)i(x,t) as small perturbations on a background flow consisting of a steady vertical vortex 

5 



U(x), with corresponding surface deformation r)o(x)] u <C U, where U denotes a typical 
value of U(x), and r]i <C r] . Substituting v± = U(x) + u(x,t) and 77 = rj Q (x) + r}i(x,t) 



into ( |2.6| ) and ( |2.4| ) leads, to leading order in the small perturbations 771 and u, to 

- gV ±V o = (u ■ V ± )U = ^V ± U 2 -Ux rotU, (2.7) 

and 

U ■ V ±Vo = 0. (2.8) 

These equations allow the computation of r]o in terms of U for the background flow. The 
first order equations are 

d t u + (U ■ Vl)u = -(« ■ V ± )U - gV^, (2.9) 

dtTii + (U-V±)Tii = -hV±-u (2.10) 
-[t7 V_l u+ (u ■ V ± )?7o]- 



Taking the horizontal divergence of ( |2.9| ), we obtain 

d t V ± ■ u + gA ±Vl = -V± ■ [{U ■ V ± )u + {u ■ V±)U], (2.11) 
where Aj_ = Vj_, and rearranging the right hand side of ( [2.1 1|) gives = 1, 2) 



AVi ■ u + gA ±Vl = ^(diU^djUi), (2.12) 

where D t = d t + U ■ Vj_. Taking the difference between D t of Eqn. (|2.11|) and h times Eqn. 
( gig ) leads to 



AV - c 2 A ± ^ = -D t [r) V ± -u+{u- V ± )77o] + 2h(diUj)(djUi), (2.13) 



where c = a/^7i is the phase velocity of shallow water waves. 

We consider the case U <C c . In analogy with gas dynamics, we call M = U/c the 
Mach number. We denote a typical length scale of the vortex by a, and the wavelength and 
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frequency of shallow water waves by A and v respectively. We will assume that wavelengths 
are small compared to vortex size ||10|| : ka = (3 ^> 1 , (k = 27r/A). 

Under these assumptions, the right hand side of ( |2.13| ) will be O(M) or 0(/3 _1 ) compared 
with the left hand side. Neglecting these terms, the final equation to be solved is 

D^T)! — c 2 A±r}x — 0. (2.14) 

Note that one might be tempted to neglect (U ■ Vj_)?7i with respect to d t T]i on the grounds 
that U <C c. However, it is possible to have (U ■ V±)i]i ~ dtTji without violating this small 
Mach number assumption by considering (as we do in Section 4 below) a background velocity 
U ~ uja with length scale a and frequency scale o>, together with v to and ka 3> 1. 

The equation ( p. 14 ) is readily obtained, under the same assumptions, in the diffusion of 



acoustic waves by a vortex ||. The physics is the same since acoustic waves and shallow 
water waves are both nondispersive, and the results of Sections [TV] and [V] are valid for both 
types of waves. They depend only on two parameters, the dimensionless wave number (3 
and the Mach number M, and when those parameters are the same the results that hold 
for the surface elevation rjx may be transposed, quantitatively and with no change, for the 
scattered acoustic pressure. 

III. ANALOGY WITH THE AHARONOV-BOHM EFFECT 



The wave equation ( 2.14j ) posesses a close analogy with the quantum mechanical wave 



equation describing the Aharonov-Bohm effect, in which a magnetic vector potential in- 
fluences the dynamics of a charged particle in a region where the magnetic field vanishes. 
This cannot happen in classical electrodynamics |§. In its simplest form, this effect occurs 
when a beam of particles with charge q and mass m is incident normally on a long thin 
cylinder containing a magnetic field B( x) parallel to its axis. The Schrodinger equation in 
the presence of a magnetic vector potential A is 

_L HW - qA{x)fi,{x) = ^4>(x), (3.1) 



where % is Planck's constant. Outside the cylinder, A(x) = ($/27rr)#, with $ the magnetic 
flux contained within the cylinder, and 9 an azimuthal unit vector. Of course, B = outside 
the cylinder. 



Both Equations (|2.14j ) and (|3.1|) allow for a solution of the form 

exp[— i(k ■ x + a9)], 

where a = vT/(2itc 2 ) = kT /2irc in the fluid mechanics case and a = —q^/2nh in the 
quantum mechanics case. This is an exact statement in the latter case, while in the water 
wave case it is approximate, because ( |2.14j ) is valid only when M <C 1 and /? > 1. Except 
for integer values of a, this is a multivalued solution. Berry et. al. |T| showed how fixing 
this multivaluedness leads to a solution that is a superposition of dislocated wavefronts 
and scattered waves. This was achieved by solving the Schrodinger equation Q3.1|) with 
impenetrable boundary conditions: ip = at the surface of the cylinder. The appropriate 
boundary conditions in the fluids case are continuity of velocity and of surface elevation. 
We now turn our attention to solving Eqn. ( 2.14J) under these conditions. One important 



physical difference between the classical and quantum mechanical cases is that in the latter 
the phase of the waves cannot be measured, while in the classical case it can. Table 1 
compares these two cases. 



IV. SCATTERING OF DISLOCATED WAVES BY A VORTEX 



As an example, we consider a scattering problem by a circular uniform vortex with 
vorticity u and radius a surrounded by an irrotational flow. Using polar coordinates (r,9), 
the background flow is given by JXTJ] 



U = < 



if r < a 



^-9 if r > a 



(4.1) 
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where T = irua 2 is the circulation. Eqn. ( 2.14 ) will be solved separately for r < a and r > a, 
and the results matched with a continuity condition. 



Inside the vortex we have, from (2.14[), 



[(d t + (u/2)d e y - c\di + (l/r)d r + (l/r 2 )d 2 e )] m = 0. 



(4.2) 



We look for solutions that evolve harmonically (with a single global frequency v) in time, 
and Fourier decompose them in the polar angle 9: 



i^nd—vt)] 



?7i = Re[^r)i n e 



where Re stands for the real part. Introducing this expression into (fO) we obtain 



(4.3) 



dr 2 r dr r 2 ^ j Vi 



0. 



| v — nu j 2 1 



(4.4) 



Equation fl4.4|) has both Bessel and Neumann functions as solutions if k n ^ 0. Regularity 
at the origin will exclude the latter. If 2v/u is an integer, k n (= k\l — n/n^) vanishes for 
n = rid = In this case, ( ^4l) can be solved by assuming r)i n oc r p . Substituting this 

into ( |4.4|) , we have p = ±n and negative values of p are excluded, again because of regularity 
at the origin. Thus we have 



rn(r,6,t) =Re 



E 



J\n\{knf) i(n9-vt) 



% & J\n\(k n a) 



+ 



(4.5) 



C{n d )a nd I- J 



r\ 



rid 



t i(n d 6-ut) 



where the a n are as yet undetermined coefficients and C(n^) = 1 when Ivjuj is an integer 
and vanishes otherwise. 

Outside the vortex, r > a, the assumption U 2 /c 2 1 reduces ( |2.14j ) to 



d l t + — d e d t - c\d 2 + (l/r)d r + (l/r 2 )d 2 ) 



Vi =0. 



(4.6) 



Inserting the form ( |4.3|) of r/i into this equation gives 



+ k 2 fj ln = 0, 



' d 2 Id n 2 + 2na 
dr 2 r dr r 2 



k= V -. 

c 



(4.7) 



with a = vY /2ttc 2 . We wish this parameter to be of order 1. Following Berry et al. [|l|] we 
write the surface elevation outside the vortex in the form 



r)i = Re(r)AB+VR), (4-8) 

where 

Vab = E K^^e^-^, m ee Vn* + 2na, (4.9) 



with (3 = ka, and 



The coefficients a n , b n and c n are defined so that they denote the amplitude of the wave com- 
ponents at the vortex boundary r = a. In order to obtain these coefficients, the continuity 
of 7] and Vj_?7 at r = a is required. This gives two relations: 

a n = b n + c n , (4.11) 

aj, .$m = kh JJM + Cn ^M. (4 , 2) 

J\ n \(k n a) \ J m ([3) H^(/3)J 
If n = rid, the corresponding relations are 

a n = b n + c n , (4.13) 

a \ J m {l3) H^((3)J 

The third condition comes from the boundary condition of rj at infinity. We require 
that the asymptotics of tjab coincides with the dislocated wave incident from the right plus 
outgoing waves only. This leads to (see Appendix A) 



J m (J3) K 

Using the notation 
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i) m (4.15) 



In = 



1 - 



na 



(4.16) 



and the relationship zZ' v {z) = zZ v _\(z) — uZ v (z), where Z v (z) is any one of the Bessel 
functions, the following expressions for a n and c n are obtained, when 7 n 7^ : 



-l) m Jm(P) 



H 1 ^) , J m _!(/3) 



+ 



him MP) 



(4.17) 



where 



-l) m Jm{f3) 



l n J\n\-\{Pln) , Jm-l(P) 1, 
! — ! _L — — ( rji 

Win) MP) P { 



H) 



^l(^) + 7n./ W -l(^n) + l (TO _ |n|)> m = v ^2^. 



M(Pln) P 

If 7 n = 0, i.e., n = rid , these formulae are replaced by 

„ _ H) md M(P) 



*n d 



K d -i(P) M-M 
HUP) M(P) J ' 



(4.18) 



(4.19) 



(4.20) 



{-i^MiP) 



{n d + m d J/5 H — -T7- 

Jm d \P) 



where 



■^§tjjr + (rn d + n d )p-\ m d = ^ n d 2 + 2n d a. 



>m d {P) 

These expressions are the main algebraic result of this paper. 
The limit r — > oo gives the surface elevation as 



i(—kr cos 0+aQ—vi) 

■ e i(kr-ut) gin ^ 



(27riA;r) 1 /2 C os(#/2) 

„i(fcr— i/t) 



-l)[ Q le 



»([a]+l/2)0 



(4.21) 



(4.22) 



(4.23) 



(27rifcr) 1 / 2 

where the function G is defined in Appendix A, and [a] denotes the integral part of a. The 
second term in the right hand side of the equation diverges for 6 — > 7r. This is because 
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this asymptotics is valid everywhere except in a narrow sector centered around the forward 
direction, 9 = n, of angular width 0(1/ Vkr), where t]ab cannot be separated into incident 
and scattered waves |I2| , and it does not make sense to speak of a forward scattering 
amplitude. This peculiarity was already pointed out by Aharonov & Bohm || in the case 
of scattering by a point vortex. 
Also 

m - ( -4-) 1/2 e« kr ~^ Etttt^^" 6 -™/ 2 ). (4.24) 



The sum of the last term of (|4.23 ) and ( |4.24|) is the correction to the Aharonov-Bohm 



scattering amplitude that comes from the matching of the surface elevation and of its gradient 
inside and outside the vortex core. 

Berry et al. have calculated a correction for different boundary conditions. They consider 
the finite radius of the scattering solenoidal field, which is considered as impenetrable. In the 
quantum mechanical context, the scattering is due to the magnetic field inside the solenoid, 
and in an hydromechanical context it could be a solid body rotating in a perfect fluid. Their 
result reads 



_^ \ 1//2 i(kr-ut) J\n-a\(P) i( n e-Tr\n-a\) 



Since the usual scattering cross section is not defined in the forward direction, it is 
interesting instead to compare the difference in the far-field correction to the Aharonov- 
Bohm wave function (obtained in the limit of zero cylinder thickness) calculated by Berry 
et. al. on the basis of Schrodinger equation, and our own calculations obtained on the basis 
of the fluids equations. The general asymptotic form of the scattered wave rj$ is 

Vs ~ f{9)r~ 1/2 e i{kr - ut) , (4.26) 

with a scattering amplitude f(9) given by 

V 27TIK 
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In the following section, we will compare the correction to the Aharonov-Bohm scattering 
amplitude for a vortex, that is 



f{6) = G(9, -tt/2) + 2 £ -^e me Hr, (4.28) 



with the correction for an impenetrable solenoidal field, 



/ Berry (tf) = 2^ l^ e^l-D, (4.29) 

n H \n-a\\y) 



V. NUMERICAL EXAMPLES 

The solutions we have obtained are parametrized by two dimensionless numbers: a = 
vY '/2ttc 2 and (3 = ka. That is, for a given incident wave, they depend on vortex radius and 
circulation as independent parameters. The Mach number is related to a and (3 through 
a = Mf3. Scaling radial distance with the vortex radius, r' = r/a, the analytical expression 
of the surface displacement is summarized as follows: 

r)i = Re rjc, < r < 1 

Vc = E a ^f^ S^ (5.1) 
+C{n d )a n / n *S n * e - vt \ 

rji = Re(r] AB + i] R ), r' > 1 

Vab = ^2(-i) m J m {Pr')e i{n6 - ut \ (5.2) 

n 

H m{P r ') i(nO-tst) /k o\ 

m ~V n Huv) e • (5 - 3) 

where m = \fn 2 + 2na. 

We have numerically computed the total surface displacement given by (|5.2[ - |5.3|) for 
several values of the parameters a and (3. In order to approximate the series in ( |5.2[ - [5.3D 
by a finite sum, it is necessary to estimate their convergence. This is done in Appendix B, 
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where it is shown that r) c is an absolutely and uniformly convergent series, and that 7]ab and 
T)r are both absolutely and simply convergent series. As an illustration, absolute values of 
the coefficients a n and Cyi £1X6 plotted in Fig. 1. 

Since convergence of the series expansions for tjab and rjn is not uniform, the number of 
terms to keep in the infinite series depends on the value of r'. In practice, we compute the 
patterns of the surface displacement in the region \x'\, \y'\ < 5[(x',y') = (r' cos 6,r' sin 6)\ 
by the finite sum of ( |5.2|) and (|5.3| ) with \n\ < 50 for /5 = 10 and \n\ < 30 for (5 = 5, but 



we keep more terms, \n\ < 90 in (|5.2| ). Fig. 2 shows the resulting displacements for (3 = 5 
and a = 0.5, 1, 1.5, 2, and Fig. 3 for /3 = 10 and the same values of a. The dislocation of 
the incident wavefronts by an amount equal to a is clearly visible. The outward travelling 
scattered wave is also visible. Note the strong interference patterns between scattered and 
incident wave. 

Another illustration is given in Fig. 4, where we substract to the total field the dislocated 
wave. The scattered wave appears as an outgoing cylindrical wave, with a clearly visible 
dislocation in the forward direction. This is the part of the wave that does not decrease as 
1 / \J~kr, and that ensures single- valuedness of the total field. Note that the representation is 
for an half- integer value of a, and the scattering amplitude is exactly zero in the direction 
= 71 0! the comparison between the two figures clearly shows the exact compensation of 
the dislocation in this direction, because of destructive interference, to yield a single valued 
total wave field. 

Finally, Fig. 5 shows the absolute value of the correction to the Aharonov-Bohm scat- 
tering amplitude, compared to the correction calculated by Berry et al. For a > 0.5, the 
parameter m is imaginary for small negative n. This induces very different partial ampli- 
tudes for exp(— inB) and exp(ind) when n is small. Our calculations thus predict a forward 
scattering with a strong asymmetry, which increases with a as shown in Figs. 5 (c) and 5 
(d). This asymmetry effect is observed both in experiments on water wave scattering by a 



vortex |n| and in direct numerical simulations of sound scattering by a vortex |TJ[ . As can 



be seen from the dashed curves in Fig. 5, this asymmetry is absent in the calculation of 
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Berry et al. For a < 0.5, the parameter m is real for all n and the scattering in the forward 
direction is much less asymmetric (Figs. 5 (a) and 5 (b)). 
All calculations were performed using Mathematica ||15|| . 

VI. CONCLUDING REMARKS 

We have computed the surface displacement due to a surface wave interacting with a 
vertical vortex in shallow water when the vortex core performs solid body rotation, the wave- 
length is small compared to the vortex core radius and the particle velocities associated with 
the wave are small compared with the particle velocities associated with the vortex. When 
the parameter a = vF /2nc 2 is of order one or bigger, the wavefronts become dislocated. The 
scattered waves interact strongly with the dislocated wavefronts and produce interference 
patterns. The differential scattering cross section is strongly peaked along a direction at an 
angle with respect to the incident direction. This is in contrast with previous calculations 
of Berry et al. |IJ in the case of quantum mechanical scattering by an impenetrable cylinder 
of finite radius. In the sequel to this paper 0, we will show that these properties roughly 
persist when the depth of the water increases. This is important because deep water waves 
are much more amenable to actual experiments. 
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APPENDIX A: ASYMPTOTICS 

In this Appendix, we study the asymptotic behavior of the function ( |4.9| ). To this end, 
we use the computations in the Appendix of Berry et. al.'s paper In order to avoid 
confusion, we use the following notations: Our definition of m is called m new = \/n 2 + 2na, 



15 



whereas the function used in fl|] is called m \d = \n + a\. Similarly, we note b n the constants 
in the series representation of tjab in our work, whereas the constants for r^g J are noted 
K- 

Solutions to our Eqn. (fOj) are Bessel functions of order m new . Moreover, the dislocated 
wave 

exp[-i(k ■ x + a6 + vt)] (Al) 



is a solution of Eqn. (|2.14j) asymptotically, that is for kr ^> a. Consequently, it is appro- 



priate to take as a boundary condition at large distances from the vortex that the solution 
should approach this dislocated wave. 
Let us consider 

VAB = 2^ b n~f 7^V^ 

n J mnevt\t J ) 

Coefficients b n should be determined from the boundary condition that t]ab should tend 
asymptotically to ( |A1|) plus purely outgoing cylindrical waves. The representation 

J m (z) = — / e l{mt - zsmt) dt (A2) 

is still valid for m = m new , even for those m's that are purely imaginary (Ref. |Tf|] p. 954, 
formula 8.412.6). This happens when a is bigger than 0.5, and for those n's satisfying 

-2a < n < 0. 

Next, we note that as n grows, with a ~ 0(1), the difference between m id and m ncw 
decreases rapidly. Consequently, there will be an N, such that, if n > N, the difference 
between the two m's will be smaller than any preasigned value. Let us write 

Vab-Vab + 2. [ b nj (8)- hn l (3) e + 

|n|<iV V u m ncv ,\y) u m \ A \F) / 

I Jm ncw (kr) ^ J mo Jkr) \ ^ in9 (A3) 



\n\>N V 



• v— 

=Rj\ 
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The wave rf^^ is the original result of Aharonov & Bohm || , and represents the scattering 
by a point vortex, hence the notation. The decomposition (|A3|) is interesting only if the last 



sum, Rn, is small when N is sufficiently large. We will see that it is indeed the case, and 
we temporarily drop it from the calculations. 
We know that if 



then rj^g 1 gives a dislocated wavefront plus an outgoing cylindrical wave. Next, if 

b n 



we may write, 



/7T-T-2O0 
dte- iM G{e,t) 
-TT+ioO 



where 



G(9 t) = ^ e in6 (j' m ^i t ~' K / 2 ) — e imold(t-7r/2)~j (A4) 
\n\<N 

is an analytic function of t, because it is a finite sum of analytic functions (exponentials). 
Also, it is dominated by the contribution from low n's. For kr — > oo, tjab can still be 
evaluated using steepest descent. Since G does not have any poles, the pole contribution to 
T]ab is the same as that of Berry et. al., namely Eqn. (A4) of ffl. This is good, since it 
is just the dislocated incident wave. Also, G(t = 7r/2) = for all 9, including the forward 
and backwards directions. This means that there are no further contributions from the 
t = 7r/2 saddle point. This is also good, since the outgoing character of the scattered wave 
is preserved. On the other hand, 

G(9, -tt/2) = £ e ine (e" im — * - e" im ° ld7r ) . 

|n|<V 

This is different from zero for all 9, including the forward and backwards directions and we 
have, outside a small angular sector around the forward direction, the following asymptotic 
behaviour at large distances: 
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p ikr 

VAB (r - oo) = rf-\r -> oo) + -_= G(0, -tt/2) (A5) 

V 27nkr 

This result differs from that obtained by Berry et. al. on the basis of Schrodinger's 
equation. 

Let us turn back to the behavior of Rn at large N. We consider the behavior of 

\Rn\ < E \(-i) m -- mM J mnc M) - Jm old (z)\ (A6) 

\n\>N 

where z = kr is a fixed number. Using the asymptotic expressions of Bessel functions for 
large values of the index (Ref. |16|, formula 8.452.1), we have 

g|n+a| (tanh Si— Si) 

J\ n+a \(z) ~ , |n + a| = z cosh^, (A7) 

'2w\n + al tanh Si 



g \/'ra 2 +2na(tanh £2 — <52 ) 

Jj „2 + o n J z) ~ Vn 2 + 2na = zcosh5 2 , (A8) 

y 2 r K\/n 2 + 2na tanh £ 2 

where ~ means that we consider only the dominant behavior at large n. An important 
point is that these expressions suppose that n > z. The following study concerns simple 
convergence of the series Rn, for a fixed value of z, not uniform convergence valid for all z. 
We define 



a 2 



2z\n + a 



0(l/n), 



so that large n behavior means small e. It is easy to show that 5 2 ~ Si — e/ sinh <5i, and that 
y/n 2 + 2na ~ In + al — ez. We then deduce that 



7T 



Vn 2 + 2na(tanh<5 2 — S2) ~ |n + a|(tanh5i — <5i) + ez^i, 

and that 



Vn 2 + 2na tanh 5 2 ~ |n + a| tanhJi — e^cosh^/ sinh^i, 



5i = Argcosh = O(logn), 

z 

so that 

ez5i = O (log n/n) <C 1 

for large n. We have thus 

(-f™" mdd fe(^) - J\n + a\(z) ~ + W2)J| n+a |(^) = O 

Using only the very rough inequality logn/n^Jn < 1, we can now conclude on the asymptotic 
behavior of at large N. Begining with (|7£6|), we obtain 

oo 

\R N \ <J2N~ n < N- N (A9) 

N 

up to prefactors that we have dropped. The important point is that, indeed, Rn is a very 
small correction at large N, which validates the preceeding analysis. As a last remark, we 
insist on the fact that all calculations are done for a fixed value of z, and that N, at a 
prescribed accuracy, may depend on z. 



logn n 



n 



n 



APPENDIX B: CONVERGENCE 



In this appendix, we discuss the convergence of the numerical series (|5.2| - |5.3| 



The simplest case is that of t]ab- In this case, the variable r' may extend toward infinity, 
and we fix its value in the calculations. Therefore we can conclude only on simple convergence 
of the series, not uniform convergence. For large n, m ~ n. Using the formula (|A7|), and for 
a fixed value of z' = j3r', the angle 6 = 0(logn) and we get 

J m (z')-^r{-)\ 

'n \nj 



so that most clearly the series ( |5.2| ) is absolutely simply convergent. 

In the coefficients a n and c n , some functions depends on j n P, and from ( |4.1(j| ) we get 



7 n /3 ~ nM where M <C 1 is the Mach number. Thus 7„/5 <^ n, so that to get the asymptotic 
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behavior at large n of J n {lnl3) we use the same formula ( ^Tf) as before, but the angle 5 is now 
a constant of order one. We then deduce the asymptotic behavior of A n from its expression 



(EH) 



A, 



^^ + 7.J W -x(^n) + l (ro _ |T||) = Q(w)- 



(Bl) 



=0(1) =0(n) =0(l/n) 

We have seen in the preceeding paragraph that the convergence of J m (/3) is extremely fast, 
which ensures convergence of a n . For < r' < 1, the term J| n |(7 ri /3r / )/J[ n |(7 n ,/3) takes the 
maximum value at r' = 1 for sufficiently large values of n. Then the absolute convergence 
of the sum (|5.2|) is guaranteed by the absolute convergence of the coefficients a n . Since the 
support of rj c is compact, this convergence is uniform. 

Rather easily, we get that the asymptotic behavior of c n is that of J m (/3), which converges 
very rapidly. Let us introduce 

j3r' = raj cosh5 1; (3 = m/ cosh5 2 - 

We have that #i < 62, both being asymptotically of order logn. Using one more time the 



formula (|A7|) , and its equivalent for Neumann functions (Ref. |T6[], formula 8.452.2), we get 

Hl(f3r') E 1 -iF 1 



Him 



E, - iF, 



exp[m(<5i — 62 



where {i = 1, 2) 



exp(mtanh5j — m5i) 

\ / 27rmta.nh.5i 
exp(m5j — mtanh^j) 



Ei = 

Fi = 

Jnmta.iah.Si/2 

which converges exponentially fast because 5± — 62 < 0. We deduce that t]r is an absolutely 
converging series. However, in this case, r' takes values in an infinite interval so that the 
convergence is only simple. 
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FIGURES 

FIG. 1. Plot of the absolute value of the coefficients a n (a) and c n (b) versus n in a log- linear 
scale for (a,/?) = (0.5,10), denoted by dots, (a, (3) = (1.5,10), denoted by empty circles, 
(a, (3) = (0.5,5), denoted by filled circles and (a, (3) = (1.5,5), denoted by empty squares. Note 
the asymmetry with respect to n — > — n 

FIG. 2. Density plot of the surface elevation for the total wave patterns for (3 = 5, a = 0.5 : (a), 
a = 1 : (b), a = 1.5 : (c), a = 2 : (d). The greyscale is linear with surface amplitude (arbitrary 
units). The dark circle indicates the vortex location. Vortex rotation is counterclockwise. The box 
size is 10 x 10 in units of the vortex radius a. The incident wave comes from the right edge of 
the box. Note the dislocated wave, and the asymmetric scattering that occurs practically within a 
single quadrant. 

FIG. 3. Same as figure 3, for (3 = 10, a = 0.5 : (a), a = 1 : (b), a = 1.5 : (c), a = 2 : (d). 

FIG. 4. Density plot of the surface elevation for a dislocated incident wave, with parameters 
a = 1.5, (3 = 5(resp. (3 = 10) : (a)[resp (c)], and for the difference between the total wave field 
with the same parameters, represented in Fig 2 (c) [resp. Fig 3 (c)] and the dislocated wave : 
(b)[resp. (d)]. Figures (b) and (d) correspond to the scattered wave generated by an incident 
dislocated wave. Such a scattered wave is itself dislocated in the forward direction, thus ensuring 
single valuedness. 

FIG. 5. Polar plot of the absolute value of the correction to the Aharonov-Bohm (i.e. point) 
scattering amplitude, in the case of an impenetrable cylinder (dashed line) and in the case of a 
vortex (solid line), for (3 = 10 and a = 0.25 : (a), a = 0.5 : (b), a = 1 : (c), a = 1.5 : (d). 
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TABLES 



TABLE I. Aharonov-Bohm effect in quantum and classical mechanics compared and contrasted. 



Quantum mechanics 


Fluid mechanics 


magnetic held 


vorticity 


B = V x A 


u> = V x U 


vector potential A 


velocity U 


magnetic flux <I> 


velocity circulation T 


wave function ip 


£ 1 " 1 j_ 

surface displacement rj 


r\'\£i\r\r'ni"\r\Y\ t~» q t pi tyi p^ot 


niQinrnf inn t~*pi vtiTYi (^i'f^v 


a = -q^/2-Kh 


a = kT /2irc 


dislocated wave is 


dislocated wave is 


an exact solution 


approximate solution 


phase is 


phase is 


not measurable 


measurable 
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